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Abstract 

We study the lowest energy i? of a relativistic system of identical bosons bound 
by harmonic-oscillator pair potentials in three spatial dimensions. In natural units 
h = c = 1 the system has the semirelativistic "spinless-Salpeter" Hamiltonian 

N TV 

H = ^\jm'^ +pi+ ^ 7|r, -r,f, 7 > 0. 

i=l j>i=l 

We derive the following energy bounds: 



E(N) = min 

r>0 



where P = 1.376 yields a lower bound and P = 3/2 yields an upper bound 
for all > 2. A sharper lower bound is given by the function where 
fx = m{N/ {■y{N — 1)^))^, which makes the formula for E{2) exact: with this 
choice of P, the bounds coincide for all N > 2 in the Schrodinger limit m ^ cxd . 

PACS: 03.65. Ge, 03.65. Pm, ll.lO.St 
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I. Introduction and Main Result 

Many-body problems form essential links between quantum-theoretical models 
and real nuclear, atomic, or macroscopic systems. However, even for nonrelativistic 
quantum theory, there are very few many-body problems that have explicit analytic 
solutions; the harmonic oscillator and the attractive delta interaction are well-known 
exceptions. In relativistic quantum theories the situation is even worse, in spite of the 
fact that the phenomenon of particle creation allowed by quantum field theory would 
suggest that there is no such thing as a one-body problem in that theory. Therefore, 
it is of considerable interest to study model A?" -body systems within the framework 
of the semirelativistic "spinless-Salpeter" equation. For this problem there exists a 
well-defined nonrelativistic limit which yields a useful consistency check. Specifically, 
we investigate in this paper the relative energy £^ of a system of N identical bosons 
represented by a semirelativistic "spinless-Salpeter" Hamiltonian [1,2] of the form 

N N 

H = J2\/^'+Pi+ E 7|n-r-,p, (1.1) 

i=l j>i=l 

where m is the boson mass, and 7 > is a coupling parameter, and we have chosen 
units in which h = c = 1. The operators are defined [3,4] in the momentum- 
space representation where they become multiplicative operators (c- variables). The 
present work is an extension to the case of N bosons of our earlier study [5] in which 
we derived energy bounds for the corresponding 1-body problem. We may compare 
H with the corresponding Schrodinger A/"-body problem with Hamiltonian 

N 2 ^ 

^s = Ei^+ E 7|n-r-,|^ (1.2) 

i—l j>i=l 

Given our goal of investigating the relative (that is, binding) energies, both of 
these Hamiltonians have the unwelcome feature that they include the kinetic en- 
ergy of the center-of-mass motion. This is easy to remedy for Hs, but a correct 
form is not so immediate in the relativistic case H. The exact solution to the 
AT-body harmonic-oscillator problem is periodically "rediscovered" but has been 
known at least since 1935 when Houston [6] solved it. Later, Post [7] studied the 
non-relativistic translation-invariant problem: the exact ground-state energy Eg 
may be expressed [8] for N >2 (in three dimensions) by the simple formula 
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Thus e is exactly the bottom of the spectrum of the 1-body Hamiltonian — A+tr^. 
In this paper we shaU prove the foUowing statement. 



Theorem 1 

Bounds on the ground-state energy eigenvalue E of the semirelativistic Hamiltonian 
(1.1) are provided by the formula 



E = min 

r>0 



N>2, (1.4) 



which yields an upper bound on E when P = 3/2, and a lower bound on E 
when P = -P(/u), where = m{N / {'^{N — 1)^))^, a function that makes the 
approximation (1-4) exact in the case N — 2. The function P{m) is monotone 
increasing with m, has bounds 

1.376 < P(m) < |, (1.5) 

and has the limit 

3 

hm P(m) = -. (1.6) 

m— »oo 2 

In the large-m limit, the upper and lower bounds coalesce to the corresponding exact 
(nonrelativistic) Schrddinger energy £^nr = Eq-\- Nm. 

The paper is primarily concerned with proving Theorem 1 . The main technical 
difficulties are twofold: to keep the fundamental symmetries of translation invariance 
and boson permutation symmetry, and to find ways of "penetrating" the square-root 
operator of the Salpeter kinetic energy. Our policy is to work with Jacobi relative 
coordinates to guarantee translation invariance of the wave functions, and to accept 
the concomitant complications of permutation symmetry. We discuss the relative 
coordinates and some of their properties in Sec. II. We shall exploit the necessary 
permutation symmetry to relate the A'"-body energy to that of a scaled and reduced 
2-body problem. The exact solution of the 1-body problem is discussed in Sec. III. 
It is well known that the 1-body Salpeter problem is equivalent to a Schrodinger 
problem with Hamiltonian —A+y/m'^ + r^ [9,10] . We take the position in this paper 
that the lowest eigenvalue e(m) of this problem, which is easy to find numerically, 
is at our disposal. In Fig. 1 we exhibit graphs of the functions {e(m), P{m)}. The 
extension of these results to the 2-body problem is treated in Sec. IV. The lower 
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bound discussed in Sec. V is rendered possible by an operator property introduced in 
Sec. II that aUows us, in a sense, to remove certain annihilation operators from inside 
the square-root operator. For the A'"-body upper bound discussed in Sec. VI we use a 
Gaussian wave function and minimize the energy expectation with respect to a scale 
variable. The calculation is helped by special factoring properties of the Gaussian 
and by the use of Jensen's inequality. The bounds corresponding to P = {1.376, 1.5} 
are depicted in Fig. 2, and the convergence of the bounds P — {P(/u),3/2} with 
increasing m (where = m{N/{-f{N - l)^)) 3 ) is shown in Fig. 3, for 2 < iV < 8. 

II. Relative Coordinates 

Jacobi relative coordinates may be defined with the aid of an orthogonal matrix 
B relating the column vectors of the new [pi] and old [fj] coordinates according to 

[Pi]=B[n]. (2.1) 

The first row of B defines a center-of-mass variable with every entry 1 / ^/N, the 
second row defines a pair distance p2 = (vi — r2)/v^, and the A;th row, A; > 2, has 
the first k — 1 entries Bki = 1/ \/k{k — 1), the kth entry Bkk = — a/(/c — l)//c, 
and the remaining entries zero. We define the corresponding momentum variables as 

N = (5-^)1pJ=s[pJ. (2.2) 

These coordinates have some nice properties which we shall need. Firstly, we have 



k k 

and similarly for the momenta 

k k 

i=2 j>i=l 



-r,)2, /c = 2,3,...,iV, (2.3) 



-p^)\ k = 2,3,...,N. (2.4) 



It follows immediately that if is a translation-invariant wave function which 
is symmetric (or antisymmetric) under the permutation of the individual-particle 
indices, then it follows that 



(*,p2^r) = (*,p2^), 2<z<Ar, 



(2.5) 
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and 

7r,2*) = 7r|*) , 2<i<N. (2.6) 

These expectation symmetries might suggest that the wave function is symmetric 
under permutation of the relative coordinates; but this stronger property is not 
generally true; it is the case for Gaussian wave functions. Moreover, Gaussian boson 
wave functions of Jacobi relative coordinates uniquely [11, 12] have the further 
factoring property that 

*(P2, P3, • • • , Pn) = 0(p2)6'(p3, • • • , Pat), (2.7) 
where (f) and 9 are also Gaussian. 

III. The 1-Body Problem 

We consider the 1-body problem with Hamiltonian 

Hi = \Jm^ + p2 + ^2 _^ e{m), (3.1) 

where, for coupling 7 = 1, e(m) is the lowest eigenvalue as a function of the mass 
m. By transforming this problem into momentum space we obtain the equivalent 
problem 

^1 = -A + + ^ e(m). (3.2) 

Since this Schrodinger problem is easy to solve numerically to arbitrary accuracy, 
we shall take the position that e(m) is "known" and at our disposal. We note that 
in the large-m (nonrelativistic or Schrodinger) limit, we have 

3 

e(m) ~ eNR(m) = m + -y- (3.3) 

(2m) 2 

We now define, for a given value of m, the (lowest) "kinetic potential" [13-15] h{s) 
associated with the relativistic-kinetic-energy square-root operator ^J w? + and 
the harmonic-oscillator potential by 

h{s)= inf (V',^V), (3.4) 

11^11=1 

(i/',\/"i2+p2^) =s 

where 'ip{r) is a wave function in the domain V{H\) of Hi. That is to say, we find 
the minimum mean-value of the potential, subject to the constraint that the mean 



RELATIVISTIC A^-BOSON SYSTEMS 



page 6 



kinetic energy is held constant at the value s. It follows that the eigenvalue may 
now be recovered from h{s) by a further minimization with respect to the kinetic 
energy s. Thus we have 

e(m) = min \s + h{s)] . (3-5) 

s>m 

It may be difficult to find the kinetic potential h{s) exactly from (3.4). Instead 
we construct an eff'ective kinetic potential hesis) which, when substituted in (3.5), 
yields e(m) exactly. We do this by changing the minimization variable from s > m 
to r > according to the following equations: 



hes{s)=r^, s = Jm^+{^^]. (3.6) 



Now, by rewriting (3.5) in terms of the minimization variable r we obtain the 
defining relation for P{m) as follows: 



eim) = min 

r>0 



(3.7) 



In fact, by inverting (3.7), we find the following expression for P{rn) in terms of 
the 1-body energy e(m) : 

( 2 (e{m) + ^e2(m) + 3m2) V 
P{m) = — ^ — ^ (2e{m) - ^/^{m)TZrr^j ■ (3.8) 



The graphs of e{m)—m and P(m) are shown in Fig. 1: both e(m) and P{m) 
are monotone increasing with m; e(m) — m, however, is monotone decreasing, in 
agreement, for large m, with the Feynman-Hellmann theorem for the corresponding 
nonrelativistic case. In the (ultrarelativistic) limit m — > we have Hi — > —A + r, 
that is to say, the operator limit is the Schrodinger operator for the linear potential 
in three dimensions, with lowest energy e(0) = 2.33810741. In the (nonrelativistic) 
large-m limit we have ffi — > m — (l/2m)A + r^, that is to say, the Schrodinger 
harmonic oscillator with energy e(m) ~ m+3/-\/2m. By substituting these "outer" 
energies in (3.8), we obtain the bounds 

3 

1.376 < P(m) < -. (3.9) 
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It is clear from Eq. (3.7) that the expression for e(m), as a function of m and 
P, is monotone increasing in P. Thus, by substituting, respectively, the constants 
P — 1.376 and P = 1.5, we obtain from this formula lower and upper bounds on 
the 1-body energy e(m). These bounds agree exactly with the bounds we obtained 
earlier [5,15] for this 1-body harmonic-oscillator problem. 

For later application to the AT-body problem, we now consider a more general 
1-body problem with Hamiltonian 

H = ^^m^ + Ap2 + (3.10) 

and positive parameters {/3,7, A}. We find by elementary scaling arguments that 
the eigenvalue e{m, /3, 7A) corresponding to the operator H may be expressed in 
terms of the energy function e(m) by the explicit formula 

e(m, /3, 7A) = (/3^A) * e im ' V (3.11) 



7A 



In terms of P, we therefore have 



£(m, /3, 7A) = min 

r>0 



1 

2\ 2 



/3|m2 + A(:^) | +7r2 



(3.12) 



For each /3>0, 7>0, A>0, this formula is therefore exact when 



P = P{li), where = m (^-^\ (3.13) 

it yields a lower bound when P = 1.376, and an upper bound when P = 1.5. As 
we shall see in the next section, the 2-body energy is obtained from (3.11) or (3.12) 
by simply setting A = 1, /3 — 2. It is an extension of this reasoning that will allow 
us, in Sec. V, to obtain also the A'"-body, N > 2, lower energy bound by using 
suitable values for /3, 7, and A. 

IV. The 2-Body Problem 

For the case N = 2 we have explicitly 



H = J m'^ + pI + J m'^ + pI + ^\ri - r2\'^ . (4.1) 
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Let ip{p2) be a normalized boson wave function. Then the lowest relative eigenvalue 
of the operator H is the infimum of expectation values of the form (ipjHip). But 
the boson symmetry of ip{p2) means that the two kinetic-energy terms in {ip, Hip) 
must have the same value. Moreover, in terms of relative coordinates, the operator 
P2 may be written 

Pi = (£1^. (4.2) 

Now, the operator tti would immediately annihilate ip{p2) if it were not contained 
in the square root. We claim that, inside the expectation value, the operator tti may 
simply be removed; this may be seen as an immediate generalization of the following 
observation. 



Lemma 1 

Suppose "^{xju) = iIj{x), then 



1 - 



dx dy 



*=(l-^)'*. (4.3) 



Proof of Lemma 1 

If indicates the 2-dimensional Fourier transform and our new variables are {p, q}, 
then we find J-'{'^){p, q) ~ 'i'{p)S{q), and, by definition, the Fourier transform of 
the left-hand side of (4.3) becomes 

{l + (p- g)2) ^ i,{p)5iq) = (1 + p2) ^ i;{p)5{q). (4.4) 

By transforming back to the variables {x, y}, we obtain the right-hand side of (4.3). 

□ 

Applying the generalization of this lemma to our problem in three dimensions, 
we find, for ip = i(}{p2), 



{iP, HiP) =\i;,\ 2;/m2 + Itt^ + 2^plj ijj ] . (4.5) 

By defining the pair-distance variable r = ri — r2 = V2p2, and the corresponding 
momentum as p = — iV^. = 7r2/\/2, we may rewrite (4.5) as 

{i;, Hi;) = (V, (2y/m^ + + ^r^^ v) • (4-6) 
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By using a formal relative coordinate r, we have thus recovered the well-known [16] 
2-body result: the minimum of the right-hand side of (4.6) is the bottom of the 
spectrum of H which corresponds precisely to the energy of a 1-body problem with 
the kinetic-energy parameter P — 2. This result may also be expressed in terms of 
the 1-body energy function e(m) by means of Eq. (3.11). Thus we have explicitly for 
N = 2 

E = {4^)K (^^y^ . (4.7) 

In the next section we shall apply a similar reasoning to the A'"-body problem; 
however, for N > 2 we obtain, instead of the exact energy, a lower energy bound. 

V. The Lower Bound 

Suppose that ^{p2, Pa, • ■ ■ , Pn) is a normalized translation- invariant A'"-boson 
wave function. Boson symmetry and, in particular, formula (2.3) allow us to write 

E < (*, H^) = N (m^ + p%) ^ + (^^^ 7 2p2^*) . (5.1) 

Now, from the definition of the relative coordinates, we have 



' N - 1 

Consequently, an application of an immediate generalization of Lemma 1 allows us to 
"remove" the operator tti from the square root of the kinetic-energy term and write 

E<N(^^,(^m^ + ^^Tri,) ' j + (^^) 7 2p^*) • (5.3) 

Adapting the argument presented in Sec. IV for the 2-body case N = 2, we define 
a relative coordinate r = \/2pN^ and the corresponding momentum p = t:n/\^. 
The expression for the upper bound to the lowest A'"-boson energy E then becomes 

E<n(^^, (m' + ^^^^P^) ' * j + (^) 7 r^*) • (5-4) 

The inequality (rather than an equality) in (5.4) comes only from the choice of wave 
function. If we find the infimum of such expressions over all normalized translation- 
invariant AT-boson wave functions, we would obtain the exact energy E; if we find 
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this minimum but without the constraint of boson symmetry, then the right-hand 
side of (5.4) wiU in general faU below E but will in any case be bounded from below 
by the bottom of the spectrum of the 1-body semirelativistic Salpeter Hamiltonian 



But this latter problem corresponds precisely to Eq. (3.7) if we make the parameter 
substitutions 

, = A=?(^, (5.6) 



N 

Thus, in view of the P representation (3.12), it is clear by choosing P = P{p), 
where /U = m{N/ {'fiN — 1)^))^ > 1.376, that we have established the lower bound 
(1.4) of Theorem 1. 

It is interesting to note that we can also substitute the A^-body values (5.6) 
for the parameters /3, 7, and A into the result (3.11) for the 1-body ground-state 
energy £(m, P, 7A) in order to obtain the following explicit expression for the lower 
bound: ^ 

E > {N\N - 1)^7) ' e (m ((77^) " ^^'^^ 

This expression — which is equivalent to the lower bound (1.4) of Theorem 1 — gives 
the exact energy and agrees with Eq. (4.7) when N = 2. Meanwhile, for all N > 2, 
in the nonrelativistic large-m (Schrodinger) limit it yields the exact A^-body energy 

1, 

^NR = A^m + 3 (^) " N^N- 1), (5.8) 
reproducing thus the old result of Houston and Post recalled in Eq. (1.3). 

VI. The Upper Bound 

For the upper bound we employ a Gaussian wave function of the form 

$(P2,P3, •••,Piv) =Cexp ^-a^p^j , a>0, (6.1) 

where C is a normalization constant. The factoring property (2.7) of this function 
and the boson-symmetry reduction leading to (5.4) allows us to write 



.2 , 2(Ar-i) 



2 




E<Ni(P,[m' + ^--^p'\ + { ^]li(P,r'(P], (6.2) 
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where the function 0(r) is given by 



0(r)=(^)%xp(-^) 



(6.3) 



Since the kinetic-energy operator is a concave function of the square of the 
momentum, we can use Jensen's inequahty [17] to move the expectation value (p^) 
inside the square root and thus estimate the mean value of this operator from above 
and write 



We shall minimize this upper bound with respect to the scale variable a > Q. 
We parametrize the basic kinetic-energy and potential-energy expectation values in 
terms of a variable r > by the following relations: 



By substituting these expressions in Eq. (6.4) and minimizing over the variable r, 
we establish the upper bound (1.4) of Theorem 1. 

VII. Summary and Conclusion 

This paper is devoted to the investigation of the ground-state eigenvalue of the 
semirelativistic ( "spinless-Salpeter" ) Hamiltonian (1.1) which governs the dynamics 
of a system of N identical bosons that experience pair interactions described by a 
harmonic-oscillator potential with coupling strength 7 . For a fixed coupling 7=1, 
we have represented the exact ground-state energy eigenvalue of the corresponding 
1-body problem, regarded as a function e(m) of the boson mass m, by a monotone 
rising function P(m), which is bounded by 1.376 < P{0) < P{m) < P{oo) = 1.5. 
Our bounds (1.4) on the energy of the A'"-body problem are expressed in terms of a 
formula which has this function P as a parameter. 

In Fig. 2 we have plotted the energy bounds corresponding to fixed lower and 
upper limiting values of P{m), namely, P = {1.376, 1.5}. In Fig. 3 we have kept the 
same upper energy bound, obtained with the help of a Gaussian trial wave function 
and corresponding to P = 1.5, but added the best lower energy bound of this type, 
using a "running" P = P{ii), II = m{N/{j{N - 1)^))3. The lower energy bound 




(6.4) 




(6.5) 
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of Fig. 3 is identical to the exact energy for the case N = 2. For higher N > 2, 
Fig. 3 shows the approach of both upper and lower bounds to the well-known exact 
nonrelativistic solution (1.3) in the large-m limit. 

A key ingredient in this analysis is the use of relative coordinates: only in such 
a framework could the upper and lower energy bounds be made to converge in the 
Schrodinger limit. This study of the semirelativistic harmonic-oscillator problem is a 
first step towards energy bounds valid for more general central pair interactions. 
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Figure 1. The monotone energy function e(m) of the 1-body problem defined by 
(3.1), and the monotone function P{m) used in our standard representation (3.8) 
for e(m); the function P(m) is bounded by P(0) = 1.376 < P{m) < P{oo) = 3/2. 
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Figure 2. Upper (full lines) and lower (dashed lines) bounds to the lowest energy 
E{m) of the A^'-boson relativistic harmonic-oscillator problem for A?" = 2, 3, . . . , 8 
obtained by employing the constant values P = 1.376 and P = 1.5, respectively, 
in Eq. (1.4) of Theorem 1. 
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Figure 3. Upper (full lines) and lower (dashed lines) bounds to the lowest energy 
E{m) of the A'"-boson relativistic harmonic-oscillator problem for A?" = 2, 3, . . . , 8 
obtained by employing the values P = P{iJ,), /x = m{N/{'j{N — 1)^))^, and 
P = 1.5, respectively, in Eq. (1.4) of Theorem 1. For AT = 2, the lower bound is 
exact. 



